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ABSTRACT 

The following aspects of continuous flow electrophoresis were studied: 
flow and teniperature fields, hydrodynamic stability, separation efficiency, 
and characteristics of ’dde-gap chambers (the SPAR apparatus) . Simplified 
mathematical models were developed so as to furnish a basis for understand- 
ing the phenomena and comparison of different chambers and operating con- 
ditions. Studies of the hydrodynamic stability disclosed that a M wide-gap" 
chamber may be particularly sensitive to axial temperature variations which 
could be due to uneven heating or cooling. The mathematical model of the 
f,eparaticn process includes effects due to the axial velocity, electro- 
osmotic cross-flow and electrophoretic migration, all including the effects 
of temperature dependent properties. 
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INTRODUCTION 

Hydrodynamics plays varied roles in the continuous flow electrophoresis 
of small particles, in some situations the suspending fluid does little more 
than carry particles through the apparatus, in others the flow is so con- 
voluted that electrophoretic separation is impossible. One of the complica- 
ting factors is the role of buoyancy forces which can destabilize the flow 
or establish an unfavorable, but steady laminar flow. To circumvent such 
problems it has been suggested that the apparatus be operated in a micro- 
gravity e-vironment where, due to the reduced size of buoyancy forces, the 
chamber could be made larger and field strength increased. Then populations 
of large biological particles could be fractionated into narrow subpopula- 
tions on the basis of unique surface characteristics which are reflected in 
the electrophoretic mobility. Such an undertaking obviously requires careful 
evaluations of many types. The purpose of this investigation is to furnish a 
basis for understanding the hydrodynamic characteristics of the chamber and 
their effects on the separation process. Particular emphasis is placed on 
the role buoyancy plays in establishing the basic flow and affecting its 
stability. 

Work began on this project in February of 1977 with the objective of 
assembling and evaluating current knowledge of the hydrodynamics of continuous 
flow electrophoresis. Four tasks were specified for the one year contract 
period: 

(1) Develop models to describe the flow and temperature fields; 

(2) Investigate the hydrodynamic stability of the flow field; 

(3) Develop a model to predict electrophoretic separation efficiency; 

(4) Review the SPAR apparatus and experiment. 
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Work on these tasks is complete insofar as it is covered by this 
contract and results are described in this report. The studies begun 
here continue under a separate NASA contract with Princeton University. 
The main part of the report is divided into two parts: SUNMARIES AND 
CONCLUSIONS, and DESCRIPTION OF RESULTS, where more detailed information 
is set forth. 
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SUNMARIES AND CONCLUSIONS 

More detailed information on the various subjects is contained in the 
DESCRIPTION OF RESULTS sections, here we simply summarize and discuss 
conclusions. 

Flow and Temperature Fields 

The temperature field enters the problem because it alters the electro- 
phoretic mobility of the particles and causes density contrasts which lead 
to buoyancy driven flows. The non-uniform temperature ^.'.eld itself derives 
from heat effects associated with the electric field and current. Although 
the field is three-dimensional, it is possible to simplify matters using 
perturbation methods. For present purposes we sought to establish the edge 

effects due to cooling through the side walls containing the electrodes 
* 

(cf . Figure 1) , the effects of temperature dependent conductivities for heat 
and electricity, and estimate time scales for thermal equilibration of the 
chamber. 

The edge effects were found to be substantial in that they extend into 
the chamber for distances of 1-2 chamber thicknesses from each side. This 
alters the mobility of particles in these regions and has a dramatic effect 
on the flow through buoyancy effects. It was also found that the effects 

of temperature dependent conductivities were substantial, the calculated 

* 

temperature rise being 701 larger for a wide-gap chamber (0.5 cm thick) than 
that calculated assuming constant properties. Although temperature relaxa- 
tion times for the fluid are only 10-15 seconds for a narrow-gap chamber 
(0.15 cm) 2-3 minutes are required to establish the steady field in a wide- 
gap chamber. 

* Representative dimensions, fluid properties, operating conditions, etc., are 
summarized on Table I, p. 72. A schematic diagram is on p. 7. 
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FIGURE 1. 


Schematic representation of an electrophoresis chamber. 
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For the flow field analytical solutions in two-dimensions were constructed 
to investigate ways buoyancy could alter the axial flow and to study edge 
effects. One- dimensional models were then developed to investigate the effects 
of temperature dependent transport properties on the axial flow and electro- 
osmotic cross flow. 

Several conclusions can be drawn from this part of the study. 

(a) It is necessary to include effects of temperature on transport 
properties. Models which ignore this, or treat matters inconsistantly, can 
be qualitatively and quantitatively misleading, especially with wide-gap 
machines like the SPAR device. With narrow-gap machines operated with modest 
field strengths (cf. Table 1) the use of 'average' values is satisfactory since 
temperature variations are usually small. 

(b) In wide-gap machines operating in a 1-g environment the steady-state 
axial velocity profile is unsatisfactory at modest field strengths insofar as 
electrophoretic separations are concerned. An earlier study by Ostrach (using 
a 'constant properties' one-dimensional model) identified a buoyancy-driven 
feature which made downflow operation unsatisfactory. Edge effects and 
alterations due to temperature depending properties accentuate the buoyancy 
feature making matters worse. Upflow, which was once suggested as a means of 
overcoming the difficulty, turns out to be only marginally better for the 
cases studied. In downflow the difficulty arises from a recirculating eddy 

in the center of the chamber; in upflow two eddys appear, attached to the 
front and rear cooling surfaces, and these restrict the area available for 
separation. A micro-gravity environment would suppress or eliminate secondary 
flows of this sort . Of course other means of eddy suppression ought not be 


ruled out. 
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(c) Experiments at General Electric using the wide-gap SPAR apparatus 
disclosed a meandering sample flow pattern thought to be evidence of the 
structure noted in (b) . Subsequent calculations made with the models de- 
veloped here showed that the actual power levels were far lower than those 
required according to the theory and thus the meandering flow must be due 
to another process. 

Hydrodynamic Stability 

In an attempt to ascertain the cause of the meandering observed in the 
General Electric experiments the stability of several chamber configurations 
was examined. Attention focussed on buoyancy driven instabilities for 
obvious reasons and investigations of other sorts of instability, e.g., 
those due to viscosity stratification or electrokinetic effects, etc., were 
deferred. Three sorts of instability were investigated: the inception of 
cellular motion due to heat generation in a quiescent layer, roll cells in 
a buoyancy driven shear flow and the effect of an axial temperature gradient 
on a fully developed flow. Critical temperature differences for the quiescent 
layer or the shear flow are much larger than those present in the experiments. 
For the vertical chamber with axial flow a new' two-dimensional instability 
was identified with an expecially low critical Rayleigh number. For con- 
ditions characteristic of the SPAR machine the critical axial gradient is 
(ca.) 0.5°C/cm. 
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Further experimental work will be required to establish whether or not 
the flow 'meandering' is a manifestation of the instability predicted by the 
current theory. If it is, then a micro- gravity environment will provide a 
means of avoiding it. Other types of instability mechanisms should also be 
investigated, however, so as to provide a comprehensive pictJ^e. 

Prediction of Electrophoretic Separator Performance 

Using the flow and temperature fields described earlier, a mathematical 
model of continuous electrophoretic separation was developed. The model (in 
brief) : 

(a) Accept, as input data the dimensions of the chamber, operating 
conditions and flowrates, transport properties of the buffer, location and 
size of the sample injection tube, mobility distribution of the sample, zeta- 
potential of the wall coating, m-mber and size of the sa le outlet streams, 
etc. 

(b) Predicts the mobility distribution in each of the sample with- 
drawal streams. 

Calculations were carried out using computer programs which have 
been tested on "model systems". Further refinements will be made under the 


current NASA contract with Princeton. 
SPAR Electrophoresis Experiment 
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Throughout the course of this investigation attention focussed on 
understanding the behavior of wide-gap machines and predicting their per- 
formance. We now have models of the flow and temperature fields and can 
estimate the elec.rophoretic separation characteristics of a given device. 
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A1 though refinements will be necessary’, the requisite 'first generation' 
models now exist for interpreting results from SPAR (or other) experiments. 
Furthermore, the effects of changes in process variables can be examined so 
as to optimize the separation. 

Final Comments 

Work begun during this program is being .rried on under a joint program 
coordinated by Dr. R.S. Snyder of MSFC. These tasks include: 

(a) Experimental studies (at MSFC) to ascertain the reason for flow 
meandering in wide-gap machines at 1-g. This will serve to prove or disprove 
the proposal that observed unsatisfactory operation is due to a buoyancy 
driven instability and assist in developing ways to circumvent the problem. 

i,o) Case-studies with the flow and separation models developed here to 
ascertain the ultimate (theoretical) capabilities of continuous flow devices. 

(c) Theoretical work to extend the capabilities of the model, to develop 
an understanding of three-dimension'- 1 effects and finite sample concentration, 
and to investigate other hydrodynamic instabilities which could limit resolu- 
tion. Establishing the limitations due to gravity and those arising from 
other phenome. . 

Frequent observations of particle agglomeration and ’clumping’ phenomena 
with cells underscore the need for an investigation of effects due to particle 
concentration. Recent theoretical studies (Batchelor, 1972) suggest substantial 
changes in sedimentation velocity at particle concentrations of a few percent 
but none of the extant studies deal with electrokinetic effects present m 
electrophoresis. 

(d) The development of experimental techniques to test the model using 
mixtures of well -characterized particles. This will include micro-gravity 
experiments where appropriate 
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DESCRIPTION OF RESULTS 
I. Flow and Temperature Fields 
Introduction 

Inside a continuous flow' electrophoresis chamber of the sort depicted on 
Figure 1 the temperature and velocity have a three-dimensional character. 

Cold buffer enters one end of the chamber and adjusts to the new' geometry 
within a distance, x g , which is given roughly by the formula (SchlichtLig, 

1960) 

x = 0.16 Re d. (1) 

e 

Here d stands for the half- thickness and Re for the Reynolds numner, u d/v : 
u Q is the mean axial velocity and v Q the kinematic viscosity. Since the 
Reynolds number lies in the range 1-5, the entrance length is relatively 
short and here the velocity field can be modelled as being fully developed 
(viz. independent of x) . As the buffer flow moves into the electrode region 
heat is added (volir’etrically) through the action of the electric field and 
the associated current, so, to limit the temperature rise, the front and 
back walls are kept cold. For reasons that will be explained later (in the 
section on stability) the adjustment iungth for the temperature field can be 
substantial. Thus, in the electrode region the three-dimensional nature of 
the temperature field alters the structure of the velocity field through its 
effects on density and viscosity. Another contributor)’ factor is the electro- 
osmotic cross-flow caused by the action of the field on the thin layer of 
charge in the fluid adjacent to the lateral boundaries. Although this 
velocity is typically much smaller than the axial velocity, it has a major 
role in altering the electrophoretic separation processes. 

ORIGINAL PAGK lb 
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The models developed here to describe the temperature and velocity fields 
take advantage of three facts: 

(i) The magnitude of the electro- osmotic velocity, w q , is small com- 
pared to u Q . 

(ii) The axial variation of the temperature is slow, 
fiii) The effect of temperature on thermal conductivity and electrical 

conductivity is approximately linear over the temperature range of 
interest: 0° - 35°C. 

These facts justify the use of perturbation methods to develop a description 
of the temperature and velocity fields. First, because of (i) , the velocity 
field can be split into two parts, an axial flow field due to forced and 
natural convection with a superimposed electro- osmotic flow. Next, due to 
the slow variation of the transport properties with axial position, (ii) , the 
velocity fields can be split into a fully-developed part f independent of x) 
with corrections added later to allow for axial structure. Finally, the 
simple linear variation, (iii) , makes the description of the temperature 
field particularly simple. 

Separate parts of the sequel are devoted to: 

A. Mathematical models for the structure of the temperature field 

a. A two-dimensional model in which the effects of temperature 
on thermal conductivity are suppressed. This provides a means 
of evaluating edge-effects due to heat transfer near the side 
wall electrodes. 

b. A one- dimensional model to evaluate effects of a temperature 
dependent thermal conductivity. 

c. A transient heat conduction model to estimate thermal 


relaxation times. 
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B. Models for the structure of the axial velocity field 

a. A two-dimensional, constant properties model provides a means 
of examining edge effects and buoyancy effects using the 
Boussinesq approximation. 

b. A one-dimensional, variable properties model is the basis for 
evaluating thermal effects and is used in the separation model 
described in Part III. 

C. Models for the electro- osmotic cross -flow velocity 

a. A one-dimensional, variable properties model provides a basis 
for evaluating thermal effects and is used in the separation 
model (Part III). 


Temperature Field 


The equation for the conservation of thermal energy is 


C u — * = V-kVT + a E 2 

P 9 x 


( 2 ) 


The symbols are: Cp - volumetric heat capacity, k - thermal conductivity, 
a - electrical conductivity, and E q - the electric field strength. E q is 
assumed to be a constant throughout the analysis. Both k and o vary with 
temperature in a linear fashion (see Figure 2) and so we write 


o = o 0 (l + o 1 0) 

(3) 

k = k 0 (l + k x 0) 

(4) 


where k Q and o Q are reference values evaluated at the wall temperature and 
0 stands for a dimensionless temperature, i.e., 
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FIGURE 2. 


Thermal conductivity and electrical conductivity of the A-l buffer. 
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T * T r + AT 0 (5) 

w 

The wall temperature is and AT is a characteristic temperature difference. 

If we transform to dimensionless variables with u as the characteristic 

o 

velocity and d the characteristic length, then 

Pe u || * v- [Cl + k 1 9)V0] + 1 + Oj0 (6) 


where Pe = C u d/k and AT = a E 2 d 2 /k . 

poo oo o 


Both kj and o 7 are less than unity: k^ is typically 0(10 2 ) while is 
0(10 _1 ), and so it is convenient to represent the temperature by means of a 
perturbation series in k^, viz. 

r* ( 1 

( 7 ) 

Substituting into (6) we generate a sequence of equations for 9^, 0^, ... 


9 = ©^ + kj0^ + k 2 ©^ 2 - 1 + ... 


Pe u 


30 


( 0 ) 


= V 2 ©^ + o.©^ + 


3x y‘0- ♦ 0^- «■ l 


( 8 ) 


Pe u jgil _ 7 2 0 (l) +ff 6 (0 + v*(0^ve^) (9) 

OX J- 

etc. 

From equation (8) we deduce that in the thermal entrance region 
3^/3* is 0(Pe *) and varies exponentially. Viscosity and density also 

depend on temperature so that a description of the temperature field to the 

order implied by equations (8) and (9) would entail an expansion for the 

axial velocity of the form 


u = u 


(o) 


+ u 


0) 


+ 


(10) 


ORIGINAL PAGEIb 

OF POOR QUALTH 


-18- 


where u^ is the constant properties solution, accounts for temperature 
variations, etc. This expansion would be used to furnish complete velocity 

fields for (8) and (9) . Because of the complexity of the problem it has not 
been practical here to attempt a solution which includes axial variations. 
Instead we have suppressed the axial structure and developed a 'zero-order' 
approximation with which we can assess the orders -of -magnitude of the thermal 
effects. An investigation of the details of the axial structure of the 
temperature and velocity fields is part of the work being done now under 
another NASA contract. 


Djo- D imensional , Constant Thermal Conductivity Model 


The major features of the fully developed temperature field can be found 
by solving equations (8) and (9), omitting the convective terms. The boundary 
conditions are: 

(i) isothermal side walls, 0 = 0»at y = ±1. 

(li) heat transfer through the side walls at z = ±H modelled in terms of 
a heat transfer coefficient, h, i.e., 

-ki-V-h 0 (T-T B ). 

3z 

Tg stands for the coolant temperature. For the zero-order field we have, 
in dimensionless form: 


0 = - v 2 0^ 0 '* + a 1 ©^°' ) + 1 


> (o) = 0, y = ±1, ^r = -K ° (0) 


( 11 ) 


8z 


, z = ±H 
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The Biot number, Bi, is hd/k ; H = h/d. For Bi = 0 the end walls are 

oo 

perfectly insulating and the temperature field is one -dimensional. For 
Bi ® the end walls are isothermal with ©^ = 0. Fourier transforms 
were used to solve equation (11J and the solution is 


Q^Cy.z) = f z [A n cosh X n z + BJsin ™ (1 + y) 


( 12 ) 


where 


A = — (X sinh X H - Bi cosh X H)" 1 

n ^2 n v n n n 1 

n 




n 


n 


nX n 


•>2 n 2 ir 2 

X n = T- - °1 


Figures 3-6 display representative features ot the temperature fields 
for narrow-gap and wide-gap chambers. Figures 3 and 5 are perspective views, 
Figures 4 and 6 are sections. A noteworthy feature is that the effect of 
side walls, shown on Figures 4 and 6, persists for a distance of 2-3 half 
thicknesses into the chamber at each side. For the wide-gap chamber (0.5 
cm wide) this distance is (roughly) 0.75 cm and for the narrow-gap chambei 
(0.15 cm wide) , 0.45 cm. Thus, due to the effect of temperature, particles 
within tnes^ regions will have a different electrophoretic mobility from 
those in the interior. In addition, buoyancy effects will be accentuated. 


* * * 


"FSSSSiSS 
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Perspective view of the 
operating at conditions 
Note: 2d = 0.15 cm, 2h : 


FIGURE 3. 


temperatur: ^ield for the narrow-gap chamber 
listed on Table I (uniform thermal conductivity) . 
5 cm. Eq. (12) . 
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FIGURE 4. 


Sections of the temperature field as shown in Figure 3. 
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FIGURE 5. 


Perspective view of the temperature field for the wide-gap chamber 
operating at conditions listed on Table I (iniform thermal 
conductivity). Note: 2d = 0.5 cm, 2h = 5 cm. Eq. (12). 
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One-Dinensional, Variable Thermal Conduativiti] Model 


A one -dimensional temperature field corresponds to a very wide chamber, 
H -* or one with insulated end walls at z = ±H. The zero-order field is 
given by 


e (0) (y) = i 


cos HjY 
cos N, 


- 1 


(13) 


N i • °i 


and perturbations are found from solutions to 


e (‘i . lilts' 0 ’) 2 

dy 2 1 1 dy 2 


(14) 


Using Fourier transforms the solution is found to be 


e ( 1 ) (y) = 


00 f — . 

= f £ 4 (n)sin (l+y)J 


where 


(15) 


f (n) = . 1 Hi s {(e (0) ) 2 } 

5 2 nMoj/ir 2 n 

Sn{(0^) 2 } = sine transform of (0^(y)) 2 . 

Temperature fields for narrow-gap and wide-gap chambers are shown on 
Figure 7, along with that for constant thermal and electrical conductivities. 
Figure 7 depicts matters ir. dimensionless form and it is seen that with the 
wide gap the temperature rise is 1.7 times that expected with constant 
properties at the conditions shown, for the A-l buffer this is nearly 30°C. 
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FIGURE 7. 


One -dimensional temperature field calculated so as to account 
for temperature dependent thermal conductivity and electrical 
conductivity, A - wide-gap, B - narrow-gap, C - constant 
properties (cf. Table I). Eqs. (13) and (15). 





With the narrow-gap, the temperature rise is so small, about 1.5°C, that 
variable-property effects are negligible. Here we also see (upon comparison 
with Figure 5) that the influence of variable thermal conductivity serves to 
alter the temperature rise. Accounting for the variable thermal conductivity 
lowers this maximum b> about 2°C for the wide-gap chamber at the conditions 
shown. 


One - Dimensional t Transient Response Model 


An estimate of the minimum time required to reach a steady thermal state 
can be made using a transient thermal model which ignores convection, since 
it tends to increase the equilibration time by adding colder fluid and with- 
drawing warm fluid from the region of interest. A rough estimate can be 
found from the characteristic relaxation time scale, d 2 /a Q . For the A-l 
buffer a Q = 1.39 cm 2 /s so that for a narrow-gap machine (d = 0.075 cm) the 
time scale is about 4 seconds; for the wide-gap machine the time-scale is 
45 seconds. To attain a condition near the steady-state generally requires 
2-3 'relaxation times', as shown on Figures 8 and 9. Data for the graphs 
were calculated using the solution to a transient heat conduction problem 
with heat generation, viz. 


30 _ 3^0 + 1 
3T ~ 3y 2 + 


+ o j0 


(16) 


with 0(x ,_1) = 0, 0(0, y) = 0. The solution given by Carslaw and Jaeger (1959) 
is 
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FIGURE 8. 


Transient temperature response for a narrow-gap chamber (one- 
dimensional, constant thermal conductivity) (cf. Table I). 

Eq. (17). 



FIGURE 9. 


Transient temperature response for a wide- gap chamber (one- 
dimensional, constant thermal conductivity) (cf. Table I). 
Eq. (17). 
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cos N^y 
cos Nj 


1 


+ 16 “ (-i) n exp( Q cos [2n+l) Try/2] 
o [4a 1 -(2n+l) 2 Tr 2 ][2n+l] 


(•) = [-(2n+l) 2 Tr 2 +a 1 ]x (17) 

t = a Q t/4d 2 , y = y*/d 


Obviously a more refined estimate couj.d be made by solving a two-dimensional 
with convection included, however, the characteristic time for equilibration 
would probably not change too much. Thus, in any experiment with wide-gap 
machine at least 2-3 minutes should be allowed for thermal equilibration. 

Axial Velocity Field 

Equations to describe the axial velocity are derived from the Navier- 
Stokes equations using an expansion described earlier 

u(y,z) = u^(y,z) + u^(y,z) + ... (10) 


The effects of buoyancy on the two-dimensional flow field are described in 
the first part of this section. In the second, where the effects of the 
lateral boundaries at z = ±H are ignored, exact solutions are possible 
which account fully for the effects of temperature on viscosity and buoyancy. 
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Two -Dimensional a Constant Trans-port Properties Model 

The 0(1) -velocity is described by solutions to 

0 = - K + N 2 - N 3 ©^ + V 2 u^ ( 18 ) 

where K is a constant (dimensionless) axial pressure gradient; 

N 2 = gd 3 /v£ Re and Nj * gd 3 BAT/v* Re, Re = du o /v o . The parameter 
describes the magnitude of the buoyancy effect while (-K + N 2 ) is a constant 
to be determined from the fact that the volumetric flowrate is independent 
of the temperature rise, since the velocity is scaled using the mean 
velocity, u Q . Using Fourier Sine transforms we find 

oo 

u^(y,z) = | l g s (n)sin ^ (1+y) (19) 

o 

with 

g (n) = B cosh 5?- z + C cosh X z + D 
°s n Z n n n 


The coefficients B , C , and D ar^ found from the relations 
n n n 


rnrH 


cosh + C cosh X_H + D - 0 
n 2 n n n 


N 3\ 


n 


D = '^i( !-(-!)") » ^ 

n nV |^ n 11 x n 

Finally (-K ♦ N 2 ) is found from _ue requirement that 


( 20 ) 


H 


f u^ (y.z)dyd7 = 4H 


- H -1 


( 21 ) 



FIGURE 10. 


Ax idi field for the narrow-gap chamber (cf. Table I). 

The velocity is almost indistinguishable from the fully developed 
parabola, 3(1- y 2 )/2 except near the side walls at : s ±33.3. 

Eq. (19). 
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FIGURE 11. 


Perspective view of the axial velocity field for the wide-gap 
chamber (cf. Table I). Downflow at g = 980 cm/s 2 . Eq. (19). 
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FIGURE 12. 


Central sections of the field shown on Figure 11. Note weak upflow 
in center and strong downflow along front and back cooling walls 
and near side walls at z = ±10. 
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FIGURE 13. 


Perspective view of the axial velocity field for the wide-gap chamber 
(cf. Table I). Upflow at g = 980 an/s 2 . Eq. (19). 
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FIGURE 14. 


Central sections of the field shown on Figure 13. Note downflow 
along front and rear cooling walls and near the side walls at 
z = +10. 
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Results of representative calculations are shown on Figures 10-14 
Figure 10 shows that buoyancy has relatively little effect on a narrow-gap 
chamber due to the excellent heat transfer. The maximum temperature rise 
for this case is about 1°C. End effects due to the no-slip condition at 
z = ±'.i are similar to those encountered with the temperature field. 

With the wide-gap chamber dramatic effects are present at unit gravity. 
Figure 11 is a perspective view, Figure 12 shows the velocity across two 
sections of the chamber with dcwnflow. Here buoyancy causes regions of 
reversed flow throughout a large part of the central section. With a 
finite length chamber this implies that a large reciroiiating eddy would be 
present in the central section and the chamt r would be difficult to use as 
an electrophoretic separation device. If the flow is reversed, the re- 
circulating eddy splits apart and is attached to the side walls, as is 
shown on Figures 13 and 14. Since the recirculating regions are adjacent 
to the side walls, the central section allows more-or-less free passage of 
fluid and any sample. Of course, any sample that migrated into the eddy 
structures would be difficult to recover. 


One-Dimensional Velocity Fields 


Effects due to a temperature dependent viscosity and buoyancy can be 
investigated easily using a one -dimensional model and, in fact, exact 
solutions can be obtained. Since variations are confined to the y-direction 
we have 


0 


K + + 


d du 
3y y cy 


( 22 ) 


with 


u(±l) * 0 
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The solution is 


i y, 


u(y) = - K 


. T dy l * N 2 1 7 l piy l dy 2 

y y o 

with I. determined from the fact that the velocity scale is the mean 
velocity. Thus 


( 23 ) 


1 1 


1 1 y, 


1 - - K 


( y i r i 

— dy x dy + N 2 - Jpdy 1 dy 2 dy 


(24) 


o y 


o y o 


Given expressions for the dimens ionless viscosity, and density, both scaled 
on values evaluated at 'he wall, it is a straightforward task to evaluate 
the integrals. In the calculation a two-teim expression was used for the 
temperature, viz. 


e - 


( 7 ) 


For purposes of '-.omparison we can evaluate the velocity field for 
uniform viscosity and buoyancy. The temperature field is 


e (0) (y) - j- 


cos N^y 


cos K T , 


- 1 


N ( - o. 


( 13 ) 


and the velocity works out to be 


u (0) (y) = | (1-y 2 ) - 


3N, 


N, 


1+1 1 
o 1 y o 1 n 1 


tan N 


1 


(1-y 2 ) 
cos N^y" 


^ + 7 - ? coSTJ^ 


( 25 ) 


For o i 


0 we recover an earlier result due to Ostrach (1976), 
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"■i^in t i-* 






'•'<« Uy-~ ) A— wg-4 


(o), , 3 

(y) =2 


(1-y 2 ) - 


N 3 

120 


(l-y 2 )(l-5y 2 ) 


(26) 


It is important to note that the magnitude of the buoyancy effect as 
compared to forced convection is contained in the dimensionless group N^, viz. 


JS1 

Vo 


c E 2 d 4 
o o 


Here 8 represents an average coefficient of thermal expansion, - p 1 (3p/3T). 

Thus, all other things being equal, going from a narrow-gap machine (2d = 

0.15 cm) to a wide-gap machine (2d =0.5 cm) alters the effect of gravity by 

4 

a factor of (10/3) , roughly two orders -of -magnitude. It follows that 
whereas gravity forces play minor roles with narrow-gap machines the situation 
with wide-gap machines is quite the opposite. 

For |Nj| << 1 the velocity differs very little from the familiar para- 
bolic profile characteristic of forced convection. For |Nj| >> 1 regions of 
reverse flow are present as illustrated on Figures 15 and 16. Qualitatively 
these profiles are similar to those derived from the two-dimensional model; 
quantitatively, however, they differ in the magnitude of the maximm velocity, 
which is higher here due to the effects of variable viscosity and buoyancy. 
Note regions of reversed flow in the center for downflow and adjacent to 
side walls for upflow (see insets) . 

The recirculating eddy present in the downflow configuration at these 
operating conditions renders this configuration almost useless for electro- 
phoretic separation. Even if the r low was steady, the sample stream would 
be deflected towards the wall where electro-osmotic effects would be strong. 



-51- 


FIGURE 15. 


One-dimensional velocity field calculated with allowance for 
variable viscosity and density (cf. Table I). Downflow, 
g = 980 an/s 2 . Eq. (23). 
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FIGURE 16. 


One- dimensional velocity "'eld calculated with allowance for variable 
viscosity and density (cf. Table I), Upflow, g = 980 cm/s 2 . Eq. ( 25 j 
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In the upflow configuration the eddys attached to the walls take up nearly 
80% of the cross section and severely reduce the region available for sample 
separation. 

At the same time, however, it must be recognized that the calculations 
presented here do not exhaust the set of configurations and operating con- 
ditions. The recirculations, for example, can be changed by independent 
control of the pressure gradient. Lowering the field strength has a dramatic 
effect on the buoyancy, since the temperature rise is proportional to E^, but 
would necessitate an increase in the length of the chamber, etc. There are, 
therefore, a number of options which remain to be investigated, the calcula- 
tions given here simply illustrate the hydrodynamic phenomena. 

Electro-osmotic Cross -flow Velocity Field 

The presence of a thin layer of space charge adjacent to the boundaries 

in the y-z plane (cf. Figure 1) along with the transverse electric field 

causes a well-known electro-osmotic flow (Shaw, 1969). In a parallel plate 

system open to reservoirs at z = ±H where H >> 1 the velocity profile would 

appear to be flat up to a very small distance from the wall (a few multiples) 

of the Debye thickness, < *) where a rapid transition occurs to accomodate 

.1 

the no-slip condition. For most purposes the double-layer thickness, k 
is so small that we can approximate the velocity just outside this layer 
using one of the Smoluchowski equations, 

w Q = - e E o ;/ vl (27) 

where w Q is the velocity in the z-direction, E , the dielectric constant, 

E q , the field strength and, the zeta-potential of the wall material in 
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contact with the solution in question. This apparent slip-velocity is of 
the order of a few microns per second for a potential gradient of one volt 
per centimeter. When the flow is constrained by walls at z = ±H the profile 
is forced to be (roughly) parabolic so as to acccmnodate the condition cf no 
net flow across the y-z plane. Although this velocity is too slow to affect 
the axial flow substantially the cross -flow interferes with any electro- 
phoretic separation by stretching the sample cross section. To provide an 
estimate of thermal effects and furnish a consistent representation of the 
velocity' field for use in the separations model a simplified model is used. 
In this model end effects are omitted except insofar as the w ’ is force 
the flow to turn round as sketched in Figure 17. The electro- osmotic 
velocity can be calculated from solutions to 


0 = 


iP + _L 

d£ 8y y 9y 


(28) 


with w = w at y = ±d. 

The pressure gradient arises from the need to balance viscous forces outside 
the double- layer so as to produce a field with no net flow. If we scale the 
velocity using w , and lengths using d, and write the viscosity as VJ Q y(©) to 
account for the thermal effects we obtain 


n - v -l d dw 

0 • K 2 dy y W 


which can be integrated to 

y 

( y 

w = K„ 


dy x - 1 


(29) 


(30) 
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The constant is found from the requirement that the net flow across any 
cross section in the x-y plane be zero. 

Although effects due to temperature dependent viscosity are relatively 
small and change the velocity only about 10% in the wide -gap chamber, these 
could be significant if one were modelling the separation of particles with 
small mobility' differences. 
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Electro-osmotic crossflow velocity: A - plan view showing recirculation 

caused by end walls, B - velocity profile with constant viscosity', 

(1 - 3y 2 )/2; C - velocity profile with wide-gap (cf. Table I). 
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H . HiTDRODYNAMIC STABILITY 

Introduction 

During the early s 1 ^es of development of the wide-gap machines large 
scale, irregular convective mixing was observed during experiments with dye 
tracers. One of the possible causes is the flow reversal caused by buoyancy 
(see Part I and Ostrach, 1976). At that time, the absence of quantitative 
data prevented a rest of this hypothesis. Later experiments at General 
Electric by H. Semon (1977) provided additional data and disclosed a 
persistent "wavering" of the sample stream at low power inputs; higher 
power levels caused irregilar mixing. However, the power le\els corresponded 
to maximum (centerline) temperatures only a few degrees higher than the wall 
(buffer) temperature, far below those which could produce the w-shaped pro- 
files described in Par* I. According to the analysis in Section I the 
dimensionless group gea o E 2 dVv 0 k o u o must be 0(10 2 ) or more for buoyancy to 
alter the velocity profile significantly; in the General Electric experiments 
it was 0(1). Changes in the orientation of the flow relative to gravity 
changed the allowable power levels somewhat but irregular flow persisted at 
voltage gradients necessary to give a significant electrophoretic separation. 

It was decided, therefore, to investigate the hydrodynamic stability of the 
flow. 

Several sorts of phenomena are included under the general topic of 

hydrodynamic stability: the inception of motion in an otherwise quiescent 

system, the transition from one steady laminar flow to another and transition 

from laminar to turbulent flow. To establish orders -of -magnitude, theories for 

the stability of a horizontal layer and of a shear flew were reviewed. As a 

result it appears that neither of the buoyancy mechanisms involved in these 
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two situations would be able to destabilize the flow. 

Previous work on the stability of stratified fluid layers lias centered 
on quiescent layers (cf. Chandrasekhar, 1961; Ostrach, 1964). This body of 
literature was combed to locate reports of special significance to the problem 
at hand which is distinguished by the combined processes of volumetric heat 
generation and fluid motion. Sparrow, Goldstein and Jonsson (1964) studied 
the buoyancy driven instability of a quiescent, horizontal layer bounded 
above and below by rigid walls and heated internally. Although a nonlinear 
profile does lower the critical temperature difference considerably, calcula- 
tions based on parameter values for the wide- gap chamber showed that the 
critical difference is about 10°C which is well above the values reported by 
General Electric for the horizontal (or vertical) configuration with flown 
Allowance for the effect of temperature on the rate of heat generation 
(which is not part of the Sparrow- Goldstein- Jonsson theory) may lower the 
critical temperature difference somewhat but an extension of this sort was 
not attempted, since it seemed best to understand the behavior of the 
vertical configuration which would be used in electrophoretic separation. 

Vest and Arpaci (1969) studied the stability of natural convection in 
a vertical slot, where however, since the base flow is driven by an anti- 
symmetric temperature , it is different from that in a continuous flow 
electrophoresis chamber. The stability with respect to roll wav^s, oriented 
perpendicular to the main flow was examined and a critical Grashof number of 
7880 was found for this sort of instability. In the General Electric experiment 
the Grashof number based on the maximum centerline temperature is roughly 10 
and, although the circumstances are quite dissimilar, it seems unlikely 
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that the observed meandering derives from a shear instability of the sort 
studied by Vest and Arpaci. 

If the meandering and subsequent large scale convection are the result 
of an instability, then it must be one where the critical Rayleigh number 
is small. One possibility is that the instability derives from small axial 
temperature gradients which result from uneven heating or cooling. In 
addition the thermal region near the entrance to the electrode section ex- 
tends over a region of several chamber half -thicknesses. This gradient can 
be estimated from an earlier equation describing the balance between con- 
vection, conduction and generation, 


Pe u(y) 


90 


i. , l9 c) 

ax 9y 2 1 


( 8 ) 


for temperature- independent properties. This equation has solutions which 
decay exponentially with :: and have the form 


-A? x Pe' 1 

e yy) 


The mode which has the snv.llest eigenvalue, X n , fixes the relaxation dis- 
tance, x g . An estimate of the smallest eigenvalue cai. be found from the 
problem where the variable velocity u(y) is approximated by the (constant) 
average velocity. A more exact calculation would refine this estimate but 
would not change the order-of -magnitude. It is found that the distance over 
which the temperature adjusts to the ohmic heating is approximately (Pe) (d) . 
Foi the wide-gap chamber operation at the conditions listed on Table I the 
Peclet number is about 50. Thus, since the temperature rise here is over 30 
degrees, axial gradients at the inlet (and outlet) are of the order of l°C/cm. 
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StatiZitu of a Fully Developed Flow with an Axial Temperature Grad lent 

The presence of cold fluid above warmer fluid (in the region above the 
electrodes with the downflow configuration and in the electrode "outlet" 
region with upflow) is an unstable configuration. To study the stability 
of the velocity and temperature fields they were modelled with a fully 
developed axial flow in a rectangular channel. The basic temperature field 
consists of an axial gradient of magnitude denoted by A, with lateral 
variations due to the balance between heat generation and conduction 
through the walls, i.e., 

T = T + Ax + AT0(y,z) (31) 

w 

If we scale lengths, velocities, etc., as before we obtain 


fr - - & ♦ ¥ 1 - 6iT > - ifW x • I ° * ,2u 


Gr 


where 


30 + Ra Re 
3 t Gr 

v 2 e + 1 

= gBATd 3 /v 2 

Grashof number 

= gBAd 4 /v Q a 0 

Rayleigh number 

■ U 0 d /”o 

Reynolds number 

II 

c 

o 

p 

o 

Prandtl number 


(32) 

(33) 


with u = 0 = 0 cn boundaries at y = ±1, z = +H. 
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The mathematical problem is to identify conditions under which tempera- 
ture and velocity fields other than the steady-state, symmetric forms exist. 

In particular we a r e interested in forms with an exponential time dependence, 
e 1 ^. The demarcation between stable and unstable flows is co = 0. Because 
of the linear structure of the equations we can show thr.t the imaginary part 
of u), I m (a)) , is zero so that the so-called "exchange of stabilities" principle 
is satisfied and any disturbance will grow exponentially. Thus, we simply 
look for conditions where u> = 0. The problem is decomposed into the sum of 
the steady parts u,0 and perturbations u and 0. 

This problem is very similar to one solved earlier by Ostrach (1955) , 

Here those results are extended to include two-dimensional effects and other 
disturbance planforms. In the early work the instability was identified 
through a degeneracy in a base flow, now we see that the disturbances are 
superimposed on a symmetric base flow. 

For the disturbance flow, u, and temperature, 0 , we have 

. A A 

VjjU = Ra u 

(34) 

0 Ra'^jjU 

where Vj. stands for the two-dimensional Laplace operator . 
u = 0 = 0 on the boundaries. One set of solutions will, of course, simply 
be multiples of the symmetric (with respect to the x-z and x-y planes) 
steady-state solutions. We are interested in anti -symmetric solutions, 
which represent nc change in the volumetric flow rate through the y-z plane. 

In general, the solutions to these equations can be written 



* . * 
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where 


u - U 1 * u 2 


'll u l - j2u l 


7 II u 2 ' 


- \ z u. 


X 4 = Ra 

The temperature is 


„ -l H , 

0 = Ra (u, - u^J , 


Solutions are 


u^ = sin qz(A^ sinh y^y + B^cosh y^y) 


+ cos qz(A ? sinh y^y + B^cosh y 2 y) 


Y? = q 2 + * 2 


u^ = sin qz(A^ sinh y 2 y + B^ cosh y 2 y) 

+ cos qz(A^ sinh y 2 y + B^ cosh y 2 y) 


y? « q 2 - X 2 


To satisfy the boundary conditions on the walls at z = ±H either 


(i) sin qH = 0, A 2 = B 2 = A^ = B^ = 0, q = rnr/H ; 


(ii) cos qH = 0, A^ = B^ = Aj = Bj = 0, q 


2n-l n 


The first condition, (i) , corresponds to a disturbance that is anti-symnetric 
with respect to the x-y plane with upflow on one side and downflow on the 
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other. This form preserves the volumetric flow rate. The second 

condition, (ii) , describes a flow with the requisite asymmetry if 

B, = B = 0. 

2 4 

Next, with (i), either B^ = B^ ^ 0 so that = 0 and ^ = tinr which 

gives 

X 4 = n% 4 (1 + H' 2 ) 2 (41) 




or = Aj = 0 so that B- = 0 and y^= i(2n - l)r/2 which gives 
X 4 = TT 4 [C2n - 1)2/4 + n 2 /H 2 ] 2 
With (ii) , A^ = 0 and = ±init so that 
X 4 - tt 4 [n 2 + (2n-l) ; /4H 2 ] 2 


(42) 


(42) 


The mode wnth the lowest critical Rayleigh number corresponds to th. 
velocity field 


* TlV t A A "V 

^ = sin -rj- cos -y (44) 

with the critical Rayleigh number 

Ra c = d + H' 2 ) 2 r 45) 


For the wide- gap chamber with dimensions noted on Table I, Ra c - 6.58 and 
with a narrow-gap, Ra c = 6.11. Using data for water at 10°C the critical 
temperature gradient is 0.5°C/cm for the wide-gap and 53°C/cm for the 
narrow- gap chambers. These results agTee qualitatively with experimental 
observation in that they disclose a great deal of sensitivity to axial 
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temperature gradients for die wide-gap machine. With a narrow-gap it would 
be rather difficult to achieve the axial gradients large enough to excite 
the instability. 

Although this explanation is consistent with experimental findings for 
a vertical configuration, questions still remain with regard to wilted or 
horizontal configurations. These cases have not been analyzed in detail 
but it is easy to shew (mathematically) that buoyancy effects coupled with 
an axial gradient will destr the unidirectional character of the flow. 

At present, however, neither experimental data or quantitative theoretical 
results are sufficient to assess this matter fully. 
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III. MODELING THE ELECTROPHORETIC SEPARATION IN CONTINUOUS FLOW 

ELECTROPHORESIS 

Introduction 

We turn now to the task of describing how a sample with a particular 
mobility distribution will be altered in its passage through a continuous 
flow electrophoresis chamber. To provide a first approximation we have 
chosen to base the model on one-dimensional approximations to the various 
flow and temperature fields and ignore, for the present, two-dimensional 
effects due to side walls at 2 = ±h (except insofar as they cause the 
electro-osmotic flow to recirculate) . Consequently, the temperature and 
velocity fields and the particle mobility are functions of y alone. This 
procedure is accurate as long as d/h<< 1 and regions near the side walls 
are ignored. The effect of the side walls requires much more extensive 
analysis and remains to be done. 

The model for separation is based on the fact that the velocity of a 
particle of a given mobility can be written as a superposition of an axial 
velocity, u(y), and a transverse component made up of the electro osmotic 
flow velocity, w(y) and the particle velocity due to electrophoresis, 
v m (y) . Thus for a particle of mobility, m, say, the velocity is 

iu(y) + k[w(y) + v m (y)] (46) 

and particles which enter with the sample at y o »z o will exit at y = y Q , 
z « z Q + L[w(y) + v m (y)]/u(y). If we denote the mobility distribution as 
yx,y,z), to represent the number d-" ity of particles with mobility, m, 
then 

3N 3N 

U(y) -fir- * My) * V m (y)] ^ = 0 


( 47 ) 
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describes the fact that the particles are conserved. The number density at 
a point x,y,z, is, therefore. 


N m (°»y»z - x[w(y) + v m (y)]/u(y)} 


and the problem is simply one of ''tracking". 

To predict the mobility distributions at the exit from the separator 
we suppose that there are a number of collector tubes of area (2d) (A) at the 
outlet plane. The mobility distribution in a given collector is simply 


Z.+A 


<Vi*i 


yz)dz 


(48) 


z. 

1 


N (z) 
m 


d 

N m (L,y, z ) u (y) d y 

-d 

f ~3 

u(y)dy 

-d 


(49) 


A computer program was written to implement this scheme; representative 
results are tabulated in Table II. 


Outline of the Computation Method 

The computation of the mobility distribution in the collection streams 
proceeds a> follows 

Main Input Data are: 

Cnamber dime sions (2d and Zh) 

Electrode length (t) 

Number of collection streams 
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Electric field strength 
Buffer flowrate 
Buffer temperature 

Constants in the linear equation for buffer electrical conductivity 

Constants in the linear equation for buffer thermal conductivity 

Buffer conductivities for heat and electricity at the wall temperature 

Coating mobility at 20°C 

Mobility distribution of sample at 20°C 

Location and size of sample stream. 

The program then evaluates the temperature field using Eqs. (7), (13) and 
(15) and calculates the local values of density and viscosity u^ing 
analytical formulas supplied. (Other relations can be used if necessary.) 
Then the local axial and electro-osmotic velocities are calculated using 
the appropriate equations for temperature dependent properties cited in 
Part I, Eqs. (23) and (30). Finally, particles on the edge of the sample 
are tracked to the outlet, using mobilities which reflect the local 
temperature, and the mobility distribution for each collector calculated. 
Numerical output includes 
Temperature field 
Axial velocity field 
Electro-osmotic velocity field 
Mobility distribution for each collector. 
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Results of two representative calculations, one with a narrow-gap 
chamber, the other with a wide- gap, are shown on Table II. General con- 
ditions are as shown on Table I . The sample contained two types of particles 
in equal amounts. Although both configurations show complete separation the 
wide-gap configuration separates the sample into two widely spaced collectors; 
the narrow- gap chamber barely separates the two fractions and if either 
fraction contained a distribution of mobilities there would be overlap. It 
was not possible to operate the wide-gap chamber model at 1-g without the 
recirculating flows described in Part II - so gravity had to be suppressed 
in the calculation. 

These calculations are intended to be illustrative of the sorts of 
results that can be obtained using the models derived here. Much more ex- 
tensive calculations are required to establish the differences in the 
separatory capabilities of various continuous flow devices. 



TABLE I. 


Parameters used in numerical calculations 
Fluid Properties (A-l Buffer) 


Buffer Temperature 

10 

°C 

Density 

1.0 

g/cm 3 

Viscosity 

1.33 x 10* 2 

g/an-s 

Thermal Conductivity 

5.82 x 10" 3 

watts/cm- °C 

Electrical Conductivity 

6.9 x 10" 4 

(ohm- cm) 1 

Coefficient of Expansion 

8.62 x 10" 5 

TO ' 1 

Thermal Conductivity Coefficient 

2.58 x 10' 3 

rc r 1 

Electrical Conductivity Coefficient 

3.12 x 10" 2 

to' 1 

Chamber Parameters 


Narrow -Gap 

Wide -Gap 

Electric Field Strength, v/cm 

70 

70 

* 

Gravitational Constant , cm/s 2 

980 

980 

Gap Distance (2d) , cm 

0.15 

0.5 

Width (2h), cm 

5 

5 

Length, cm 

16 

10 

Volumetric Flow, an 3 /s 

0.35 

0.7 

Average Velocity, cm/s 

0.467 

0.28 

AT r , °C 

3.17 

35.2 

Re 

2.63 

5.26 

N 2 ORIGINAL PAGE lb 

8.88 x 10 2 

1.64 x 10 4 
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n j 

2.4 x lo' 1 

4.9B x 10 1 

h 

8.2 x 10' 3 

9.1 » 10' 2 


9.9 x 10' 2 

1.1 


* 


except where noted 
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TABLE II 


Red-Blood Cell Separation in Narrow-Gap and Wide -Gap Machines 

To demonstrate the use of the separation model two runs were made at 
conditions shown on Table I with samples made up of equal amounts of 
particles ('red-blood cells') with mobilities of 2.15 pm-cm/v-s (m=l) and 
4.15 pm-cm/v-s (m=2)(at 20°C) . 

Other conditions were: 

Coating mobility 2.15 pm-cm/v-s 

Sample inlet size 0.05 cm 

Number of collectors 100 
Locations of Separated Particles 
Narrow- Gap Chamber 


Chamoer # 

1-13 

14 

15 

16 

17 

18 

19 

20 

21-100 

m = 1 

0 

.177 

.590 

.233 

0 

0 

0 

0 

0 

m = 2 

0 

0 

0 

0 

.159 

.477 

.353 

.021 

0 

* 

Wide-Gap Chamber 









Chamber # 

1-20 

21 

22 

23 

24-29 

30 

31 

32 

33-100 

m * 1 

0 

.043 

.546 

.411 

0 

0 

0 

J 

0 

m = 2 

0 

0 

0 

0 

0 

.282 

.610 

.108 

0 


* 


g = 0 
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CCMPUTER PROGRAMS 


Four programs were written in FORTRAN IV: 

ZTIMP - the two-dimensional temperature field, Eq. (12) 
VELO - the two-dimensional velocity field, Eq. (19) 
TRANST - the transient temperature field, Eq. (17) 

TEMP - the separation model. 

Listings of these are given on the following pages. 
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1J1 FCPKAT (* 3’ ,9X, 'Z' , 1JX,'lEi1?1 bATUfil') 


3313 


0006 


13C ?CiBM (' ') 


0320 

<_ 

C307 


103 ?C?",AT ( 2 (2X, 1PE 14 .b) ) 


0330 


0CC9 


15 J 5 C 5 R A I (3E 10.3,315) 


3340 


j;:9 


MAC (5, IS’') SIGBAl,EX,N t NT2BHS,NYPIS,NZPlS 


u350 

< 

■:oio 


?1 = 3 . 14 15 92t> 


03o0 



c 



>,370 



c 



0 3b0 

. 

.Oil 


DY*1./(NYPTS-1.) 


0390 


0012 


C2=«/ (NZFI3-1.) 


0400 



c 



34 10 

1 


c 

THE CUl-fi LCCF LNCSEaENTS Y , THE NEXT INNSB.10SI LOOP INCBEMtBTS 

z. 

AND TH0420 



w 

I NS . BAGS a INC LOOPS TAPE THi SERIES PCS 3 EAPEBATU«£ OJI 10 KTEfiBS 

TLB HS. J43J 



c 

THE PRCJOAI S10FES Thl lEKPEPATUBi IN THE ARBAY IFSP(I,J), Nils 

INDICES e 04 40 



c 

J XIFEfinlNG Id THI Y-LUiCTICN cdQSLINATF ANC Z-DIRECTION COOBLJNATt 3450 



c 

BFSEECUVilY . TSE ACTUAL Y ANL Z COOED! NAILS ABE SXCRiD IN AoBAYS 

IX Afc’Di 463 



c 



04 70 

l 

CO 13 


Y*-EY 


04t>0 


0014 


£C 3C 1“ 1, NYP1S 


0490 


C015 


Y*Y+LY 


O500 

i 

0016 


2*f I* ( 1 .♦ Y) /2. 


0510 

« 

CC17 


2=-CZ 


o52C 


0 J 18 


DC 20 J-1.NZPIS 


0530 


0319 


Z=iNLZ 


j5u 3 


4 i A 3 


sua*c. 


0553 


0021 


CO 10 Ml, MIS IS 


35b J 
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09/ 1 J/4D 


0022 



N*2*K-1 

0570 

0023 



LA*BD2«N*N»*I*FI/4.-5I3NA1 

05o0 

C524 



I* (1.- (-1.) **N)/N/LA»£D2 

C590 . 

0025 



LABECA*LASBC2**0.S 

0 600 

0026 



AE1*llBtCA* (2-N) 

3613 

0027 



AE2—LABBDA* (Z*N) 

0620 

0026 



Af 3*-2.*lABEDA*N 

J6 3 3 

0029 



IE (Ari*5u.) 1,1,2 

0640 

0030 


1 

1 1-0. 

0650 

C031 



JC TO 3 

0660 

0032 


2 

E 1 = £XP (All) 

367 0 

3033 


3 

IE (AE2+ 5C.) 4,4,5 

3680 

0034 


4 

E 2 s 0. 

0690 

3 035 



CC TO 6 

07 c 0 

0036 


5 

E2«tXP (AE2) 

0710 

C 037 


6 

11 (A£3*50.) 7,7,6 

0720 

0036 


7 

13*0. 

0730 

C039 



GC TC 9 

C740 

0040 


8 

13-EX? (AE3) 

0750 

0041 


3 

CONTINUE 

376 3 

0042 



I * (ZUE2)/ (LATEIJA* (1.- 2 3) - el* (1.*jt J) ) 

0770 

3043 



t<»T*Bi*T1 

37b J 

0044 


1C 

SO i-SUl + E*SIN (S*X) 

■37 ->0 

0045 



TIB? <J,I)=2.*SJ«/PI 

jflo j 

0046 


2 C 

Z2 (J)*Z 

06 10 

0047 


30 

YY U)*Y 

<jc2 J 

0 048 



lE.iil-I'.VE (1,1) 

oo30 


0 



3H40 


c 

? LCI THE TZBEiCATUEi. PI iLC AS A SOuSACi VIE lit. C IN PERSPECTIVE. 

Cb50 • 


c 



co 60 

0 049 



A* 2*NY ETS 

3870 

0 350 



CC 60 1*1, NY PTS 

uo80 

0 051 



MK-1 

■:H97 

0052 



l*2*'.z;:s 

09 '30 

0053 



SC 53 0*1 , SZPia 

3910 

C054 



L*L-1 

0920 

0055 


50 

I tip (L,K)*ltBP (NZPTS* 1-J, -.IF: 3*1-.) 

a 4 30 

0356 


6 3 

CONTINUE 

0 94 0 

0057 



SYa*iiYPTS- 1 

3953 

oose 



NZa=NZEIS-1 

0960 

0059 



K*2* NY PTS 

0970 

0960 



CC 62 1*1 , NYPT5 

3980 

0061 



K = K-1 

0990 

0362 



tC 0 2 J* 1 , NZE 

1 303 

0363 


62 

7 (tiZEIS-J,K)*TEEP (NZFTS*J,K) 

1010 

0064 



NY=.*NY?IS-1 

1 323 

C 065 



NZ**» N3PTS- 1 

1030 

0066 



CC 65 0*1 , N7 

1 340 

0 367 



CC 65 I* 1 , NIB 

1050 

0368 


65 

T.flP (J,NY?TS-I)=IiflP(J,WYPT3*I) 

1 060 

0069 



z*-rz-4 

1070 • 

0070 



CC 70 J=1,8Z 

IOoO 

JC71 




1090 

0072 


7 C 

L id ( J ) * Z 

1100 

3071 



Y — lY-1. 

1110 ' 

0074 



CC aO 1*1, NY 

1120 
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iv a itviL 

21 HAIR DATE » 78110 09/10/45 


0075 



Y*YUY 

1120 

0076 


a o 

YY (I) «Y 

1143 

0077 



JI»1F4F1* 1C. *1. 

1 150 

non 



aconrB«jT#2 

116) 

0079 



T8AX*JT/10. 

1170 


c 



1180 


c 

F 61 JIT In I FIELC FOB Y*0. 

1190 


c 



12J) 

0383 



'2 fill F (6,133) 

1210 

0081 



KBIT: (6, 110) SIG.l A 1 ,BT EBBS, B ,31 , HYPT3 , RZPTS 

1220 

CC82 



Bfilli (6,131) 

1230 

0083 



BBITI (6,120) 

1240 

0084 



tC 43 J« 1 , fii, 

1250 

0385 


43 

iiSITi (6, 10J) ZZ (J) ,TE.9P (J, NYPTS) 

1260 


c 



1273 


c 


FEINT 1H? FIELD FOB Z*C. 

1280 


c 



1290 

C086 



HfilTE (6,1)2) 

13C3 

0387 


102 

FCR8AT (• 3* ,9X, • Y*, 13X, 'T I RPifiATJ RE ' ) 

1310 

3388 



CC 41 I»1, NY 

1320 

4 339 


4 1 

Brill! (6 , 1 J J ) YY il) ,TZ8F iNEPTS,!) 

1330 

0 390 



CALL PF^VU- (11,ai,YY,lE"E,8*;.5,4.,l.'1Ai*1.5,201,131,B2,HY,C, 

1340 



X SCCSI5,J,-. 1,TaAX,u,3,C,9. f 4.,-1,17,'Ti..’1?iaAX0r£ FliU') 

1350 

0391 



SICr 

1360 

0392 



INC 

137 3 
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( 





C 



cMc 




c 



OO^J 


< 


c 



0030 




c 

Ml 


>cT40 




c 

• 


0050 


c 


c 

* 

mo ELTEFBI NcS THE DGNNBABD VELOCITY III THE ELECTBOP HOBBS IS CEL* AS A0067 




c 

* 

FUNCTION OF BOTH I AND Z. Tl'ROCGh A CALL TO THF PBOCEDUhE PEJVUE, IT 

0070 




c 

• 

THEN PLOTS A PEKSFcCTIVt VI cB OF THE VELOCI1I SOBI ALE. 

woec 


f 


c 

• 


0o90 




c 

• 

SIOSAI -PAR AHETEN I* THE EXPRESSION FOB THE ELEcTElCAL 

J1O0 




c 

* 

CCNDUCTIVITI 

0110 


f 


c 

* 

LI -I HI jIOT NUHBEB 

) 1 20 




c 

• 

U -THE ASPICT EATIO OF TH r CHANNEL 

0130 




c 

♦ 

N'.ERNS -THE NOSELF OF TERNS TO Lr TAacN IN THE SERIES FOR THE 

0140 


f 


c 

• 

VELOCITY 

0150 




c 

♦ 

NYPT3 -THE NUBHEB OF E3UALEI SPAv.cE oOBSTANI-T LUES FBOB (ABO 

0160 




c 

• 

INCLUDING) 1*0 TO 1*1 ALONG NH1CH THE VELOCITY IS 

0170 


f 


c 

• 

EVALUATED 

3183 




c 

* 

NZPTS -TH L hUBEEF OP FwUAMY SPAwED CONSTANT-! LISES FECfl (AN L 

'■•190 




c 

* 

INtiUCiNS) Z*0 TO Z^V ALCNG ■hXCH VtLOcITY IS EVALUATE^ 

c2 00 


< 


c 

• 

0210 




c 

t * 


>0220 

"4 



c 



325 3 


( 


c 



0240 



0001 



C I Hi NS ION 11(201, 1 J 1) ,YY (1 O’) ,ZZUJ1) 

J25 3 



•'orz 



SEAL LAitCA,LANBl)2,N2,l!3 

0260 


f 

0003 


100 

FCrHAI ( 1 H 1) 

327 D 



34 


111 

FCRSAT (•-', 2X , ' SIo H A 1 * ' , 2 10 . 3, 2X, • EIOT JO. *',110.3,21, 

'283 






•ASPECT f«nc * ' ,E 10. 3,2X, ' N2 - ' , c 1c. 3, 2X, * N j *',c.10.J) 

029 0 


( 

0005 


111 

ROSS AT ('0 ,2X,'NTch(!S - ' , 1 5, 2 a , * NYPTS = ' , 15 , 2 X, ' NZPTS **,15) 

cJOO 



0006 


120 

ROR.1AT ('J',9A,'Z',11i ( 'S VELOCITY*) 

0310 



0007 


130 

FCEHAT (' •) 

C 320 


( 

CCCd 


140 

FOBS AT (2 ( 2X , 1 P 1 1 4 . e ) ) 

)33 J 



1 1 '9 


150 

rCI.SAT (5 F K.3,315) 

0340 



00 10 


160 

®ORnAI ( , u , ,9X, , Y , ,11X,'X VELOCITY') 

0350 


( 

001 1 



SEAL (5, 15). 1:40=999) SJGSA1,ol,».N2,NJ,NTEB*.S.k'. TS.NZP1S 

sjf 1 

• 


0 012 



HU.i (6,100) 

0J7. 

J 


0)13 



will • (6.11)) < 5I«BA1,L1,'«,N2,NJ 

cje: 

£ 

t 

0014 



• £IT t (6,111) NTch.*E,NYP7S,NZPIS 

0390 


0)15 



OHIN-0. 

0400 



0016 



U.NAX* 0 . 

0410 


( 

0)17 



Fl*3. 1415926 

C420 



0018 



RI2=P1*PI 

0450 



0019 



r 124 - F I 2/4 . 

044 0 


( 

0 020 



f 14=FI/4. 

0450 



0021 



F I B 1 P 1 * J 

3403 



0022 



H2I2*Ei*B/2. 

0470 


€ 

0023 



7PI-2./PI 

3480 



0024 



SUH 1 * C . 

0490 



0 025 



SUNT* 0. 

.SCO 


1 

0026 



SO N 3* C . 

0510 



0027 



K=- 1 

j52 3 

> 



c 



0550 

_ 

( 


c 

Ji 

1 EB .11 NE Till CONSTANT, K, FnUfl THE 1 NT cGRAX 20 VELOCITY. 

0543 




c 



0550 



0028 



CC 5 N* 1 , NT i BBS 

0560 


t 




PS? 
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C 
f 
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r 


i 


i 


i 


< 




t 

i 


POETBAN IV G LLVil 21 MAID DATE - 76110 39/14/24 


0029 

k*m 2 

3573 

00 .0 

L A.1BC 2*A *K ’PI24-SIGR A 1 

0560 

0 33 1 

L A«£CA*1.A!1£D2*»U. 5 

3593 

0032 

?N2*K*K 

7630 

0033 

P N4* FN2*f M2 

0610 

0334 

El43f:/K 

V62C 

0335 

F7A.*n1*TANE(/*SPI2) 

0o30 

3336 

PTANH2*TANF 'l»1cuA*WJ 

0640 

0037 

S1“(21*FTAliHl-B//Flt4 

0650 

3038 

s2*5Vt*eEr2 

( 660 

003'J 

S3* ( ( P1*?TANH1-F‘XANH2/LA!!tDA)/lLAB3i;A*?TAIia2-BI))/Fll2/LAfl8D2 

0670 

3 34 0 

SUflUSUfll + SI 

766 0 

034 1 

sua2*su?u*s2 

0690 

0342 

jJMj=0013+3? 

37 0 o 

rj4i 

5 COMIkOE 

u7 10 

0044 

A1«SUM*32./i/(PI**4) 

0720 

0045 

A2*SU.12*32 ./B/ (Pl**4) 

0730 

3 346 

A3=30B3*e./SIGSA 1/FI/PI/a 

0740 

7147 

«iv=( 1 .♦i.2*Al-N3*A2-Ji3*BI»A3)/A1 

0750 

0343 

w?Ii: (0,H2) F K 

J7o 3 

0349 

112 FCSf.AT (• ‘ ,21, El v. 3) 

07*»0 


C 

3763 


C TVAILATi. VELOCITY «T IACH (Y f 2) PCI NT TAKiA. 

07*3 


C 

vsOO 

0 350 

7 Y* 1 ./ (.SYrTS-1 .) 

OslO 

305 1 

IZ-i / (NOP IS- 1 . i 

.32 7 

0352 

Y = -CY 

. 3 30 

0053 

E«1=FK-N2 

4«4 3 

3 354 

f 1*6. /PI2 

0850 

3055 

F 2= 2 . *H3*EI/5 I5H A 1 

jdb 3 

005b 

EC 30 I* 1 , AYPTS 

0s70 

3057 

Y - Y+DY 

3563 

035d 

i~il* tl . ♦ Y) /2 . 

0850 

3 359 

2--E2 

09 oO 

27e7 

10 2 r J 1 1 , NZF1S 

0910 

00 b 1 

l=i.*U 

0920 

3362 

SUS= 3 . 

o 330 

3063 

< = - 1 

09 4 0 

•J 364 

OO 13 Si*1,MLfcaS 

0950 

00o5 

K*K *2 

0960 

0066 

LAflSE2«X*X*PI24-SIGNA1 

397 3 

7067 

L*..'12CA*LA:iI>D2**7.5 

*«60 

3Co3 

7F3*K**3 

3990 

C 069 

FIANH*TA!IH (L ABBE A*H> 

1000 

0373 

Cl* FI* (rKW!«3/lAfloE2)/PK3 

1 310 

r 37 1 

C3*f2/K/LAB£D2 

1020 

0072 

C 2“C 3/ (LAflbDA*FT AK 6 - 6 I) 

1333 

0373 

AEG 1*K* (2-fc)/TPI 

1040 

0374 

A3G2*-K* (2* ■ ) /TPI 

1 353 

4 375 

ApG3*-R*FIK 

1 0 to 

03. o 

A3w4*LAflbCi* <2-« J 

l.'ij 

0377 

LOG 5--LAS6DA* („*H) 

1060 

0076 

AaG6*-2. *LA3EEA*W 

. 393 

3379 

If (ABG1.LE.-100.) E1*C. 

1 1 vO 

0 360 

IF (ARG1.UT.-133.) E 1* LI P ( A aC 1) 

1113 

0 361 

IF (A3G2.LE.-1CO.) £2*0. 

1120 
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0082 



if (ABU2.GT.-0OO.) E2*SXP (AnG2j 

1130 

J 293 



IF (ARU3.EE. -100.) EJ*J. 

114 1 

0084 



IP IABGJ.GT.-100. ) E3=£XP (AI.13) 

1150 

0085 



IF (ABG4.tE.-100.) E4*\ 

116 0 

0086 



IF (ABU4.GT.-100.) £4*£XP (A5G4) 

1170 

JC87 



IF (ARU5. tk.“ 10 0. ) E5- 3. 

118 0 

0088 



IF (A5G5.GT.-100.) E5 = EXP(AaG5* 

1190 

0089 



IF (A BU6 .LE.-IOO.) Lb- 0. 

12 J 0 

0 09C 



IF (AHGb.GT.-IOO.) £6=FX? (AHG6) 

1*10 

0 091 



F= (F1*E2)/(1.*E3) 

122 0 

C 0 92 



A 1*C1*E 

1*30 

03*3 



A2=C2»E 

1*4 0 

0 394 



A3=-C3*(£44i5)/(LA.1BDA* k 1 .- -6) -oi* (1 . *26) ) 

1250 

C 095 



A 4=-C 1 

126 0 

“2 96 



SUt x Sv)H4 (M*A2»*3*A4) *516 ih»X) 

127 0 

0097 


10 

CCNTIRUE 

1280 

0098 



0 (J.I) *-2. *SOB/PI 

1290 

0099 



IF k J (J.I) .LE.USIR) OBIR=0(J,I) 

1300 

3100 



IF (* (J, I) .G2.0BAX) uHA7. = U (J, 2) 

1 J 1 

J101 


2 C 

ZZ V J) = 2 

1320 

3172 


3 0 

n i:> = v 

1 j j 7 


c 



U40 


0 

?:c 

i, Til- V 2 I CC I" * rELl 1 A8 A SGiFAi-i VIE. EC IS HtTSP: CXI V 3. 

1 j 5 3 


C 



1 JC4 

0 103 



K=2*NYPTS 

1 j7 0 

013 4 



EC o' I-1,iiY?'I3 

libj 

0 1C5 



* = K - 1 

134 ) 

3 136 



:-2»N2fIS 

1 4 J J 

0 1o7 



Lo 5 C 0=1, 92PES 

1*10 

0 1C8 



E=L-1 

14^0 

0109 


5 0 

J (L.,K)=’J (A-.?TS*1-J,RYPTS*1-I) 

14 30 

0110 


*J ? 

CCfcl iKUl 

1 ** 4 0 

0111 



KY.\=.iYi7S-1 

• *50 

0112 



.'•:.i=n 4 .?ts- i 

j 0 

31 1 3 



K=**ST»IS 

..70 

0114 



Zo b 2 I*i,SYPT3 

1 4C J 

0115 



K -K- 1 

1*90 

Olio 



CO 62 J* 1 , RZB 

15'. 2 

0117 


62 

0 (BZPTS-J.K) *0 (IZPTS+J ,K) 

1510 

0118 



SZ*2*BZPTS-1 

1520 

0119 



RX*2»MPTS-1 

1530 

0120 



CO 65 J* 1, IZ 

1540 

0121 



CO 65 1*1, »TB 

1550 

0122 


o5 

0 (J,RTPTS-I)=U (J,HIPTS*I) 

15o 0 

0123 



U (1,1)*0. 

1570 

0 124 



9 (N2, 1) =0. 

153 0 

3125 



J (1,:<X)*C>. 

1390 

0 126 



J 4N2.RYJ-0. 

loGO 

0127 



2=-L2-W 

1610 

0128 



EC 73 J=1,SZ 

lo20 

0 129 



2-Z*CZ 

1o j0 

0100 


7; 

Zl (J) *Z 

H40 

0131 



I=-GY- 1. 

1650 

3 102 



0C 80 1*1, HI 

lbbO 

0133 



i=Y*LI 

1o7 0 

■3 '34 


3' 

YY (I) *Y 

loo J 
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0135 


UBIlc 

(6,12*) 

1690 

0136 


iiRITt 

(6,130) 

1733 

0137 


CC 4 ' 

1*1, AZ 

1710 

0138 

4 J 

asm 

( 6 , 1 4C) zz(i) ,u ( i, hipis) 

1723 

0139 


8 FII 5 

(6, 160) 

1730 

0140 


88 lit 

(6,130) 

174 3 

0141 


:c 93 

1*1, ST 

1753 

0142 

9u 

8?IIe 

(6,140) IT (I) ,0 (NZPTS, I) 

1760 

314 3 


jam* 

1 . 1 *U SIS 

1770 

0 144 


uiax* 

1. 1 *0*AX 

1763 

3145 

X 

CALL 

f£5Viii(11.2Z,TT,U,8*.5,3. , 3. , 2 3’ , 1 3 1, SZ.I I, 3, 1 3, 3, UBI N, J8AX, 
0.0,0, 9., 9.,- 1, 14,* VILCCITT F1UD‘) 

1790 

1900 

0146 

0147 

99S 

STCP 
t NC 


la 10 
1623 
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< 


f 
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c 


o £ 1 J 



c 


v „ j r 

( 


r 


► j)j) 



c 

# 

) J 4 ) 



c 

• TSANSNT CITIFSJNiS 1H £ AVERAGE TEMFEnATUcE AND CENTES T 1* f Et AT'Jit Or 

U Ji J 

f 


c 

• THE EIECT8CPH.4ESIS CHANNEL AS FUNCTIONS OF TIBI. THE fE'.JIRiD 

JOO ) 



c 

• r AIA IS II'.TEC IN ThF CUTEOI CF THE FBCGPAB. ALL THE xN.JI VAE.Ao..O 

-.7 ) 



c 

• SILF-*XP’ NATO?! IXC EFT ’*»l FCLLCNING- 

3 JH J 

c 


c 

• 

Jv. J J 



c 

• NUNTAC - TH! NUEeEB or TINES A i NhlCd TI3Pt9ATJfif iS tVAIOAIit. 

ol 30 



c 

* CTAU - TUI DIEj.NSIoNi.ESS TIBI 1C BhlCH TEflFEEATOoa 13 IVAlUA tE 

. 11 ) 

f 


V 

• NTS? NS - TH: NUflot 5 OF TE6.1S TAKEN FCP TA’J * 0 IN Thr SERIES EXh 

0 1.3 



c 

• TENFaiA lUaE. FOa TAJ G&EATEa THAN ZIBO, laE Ko.13.-H or 11. 

... 1 j' 



c 

* TAKEN IS CCNIRCLUC a! KoKLio. 

o 1 4 0 

t 


c 

♦ 

0 15' 



c 


*0 loj 



c 


317C 

r 

CCC1 

c 

tIN.KiiC.% CTE*?(1CC1) 

O 1 t V 

( 

:::2 


SIStNSICK AVIESF (1301) 

^ £ . . 

J J0 3 


lI.NtSSiSN 7 ( 1 QC 1) 

- i 1 VJ 
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